Calibrating Dark Energy 



Roland de Putter & Eric V. Linder 
Lawrence Berkeley National Laboratory & University of California, Berkeley, CA 94720, USA 

(Dated: August 1, 2008) 

Exploring the diversity of dark energy dynamics, we discover a calibration relation, a uniform 
stretching of the amplitude of the equation of state time variation with scale factor. This defines 
homogeneous families of dark energy physics. The calibration factor has a close relation to the stan- 
dard time variation parameter Wa, and we show that the new, calibrated Wa describes observables, 
i.e. distance and Hubble parameter as a function of redshift, typically to an accuracy level of 10"'^. 
We discuss implications for figures of merit for dark energy science programs. 
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I. INTRODUCTION 



Understanding the nature of the dark energy accelerat- 
ing the cosmic expansion is one of the premier questions 
in physics. The answer offers the possibility of deep in- 
sights into the nature of spacetime and gravity, extra 
dimensions, the quantum vacuum, and possibly the uni- 
fication of gravitation and quantum physics. Precision 
mapping of the expansion history provides one path to 
characterizing the dark energy, in particular its equation 
of state and time variation. 

Guidance from theory is useful to predict observ- 
able signatures for cosmological probes such as distance 
and Hubble parameter measurements, in particular what 
level of accuracy is required to distinguish between mod- 
els. From a model one can predict distance-redshift rela- 
tions etc. but the number of models is vast; one would like 
to identify model independent or at least generic charac- 
teristics of the dark energy. Indeed, such properties exist, 
as discussed in detail recently by [Ij] , for classes of behav- 
ior in the early time evolution of dark energy, valid for 
z > 2 when the dark energy does not strongly affect the 
background expansion. 

In this article we seek to extend characterization of the 
dark energy properties in terms of the equation of state 
to the entire observable history. This requires a different 
approach, calibrating the evolution through a "stretch" 
relation between the amplitude of the time variation and 
the time variable or scale factor of the expansion. The 
calibration then provides a physical basis for a compact 
and highly accurate parametrization of the dark energy 
influence on observables. 

In iJTT] we examine several diverse models, looking for 
similarities and distinctions. We introduce the cali- 
bration in mil] and discuss its relation to a standard 
parametrization of the equation of state. ijlVI examines 
the utility of the description and shows that it achieves 
robustness and accuracy at the lO^'^ level, sufficient for 
next generation data. We discuss some implications for 
figures of merit of dark energy science programs in f|Vl 
Those readers wanting to get right to the results could 
start in the middle of mill 



II. DARK ENERGY DYNAMICS 

By examining the behavior of a diversity of dark en- 
ergy models representing different physical origins, we 
can explore common and distinct elements within 1) a 
model as the parameters vary, 2) a family of models with 
some related property, and 3) different classes of mod- 
els. Families of models might consist of those with sim- 
ilar functional forms, e.g. polynomial potentials, while 
classes might be those with similar early time behaviors, 
e.g. thawing models or freezing models [2]- 

We choose five representative families ranging over dif- 
ferent physics and different evolutionary histories. These 
are the pseudo-Nambu Goldstone boson (PNGB) model, 
or cosine potential, that thaws and moves away from an 
early cosmological constant state w = —1, the family 
of polynomial potentials, also thawing, the supergravity- 
inspired SUGRA model that has early tracking behav- 
ior and then moves toward the cosmological constant 
state (freezing behavior), the modified gravity model of 
DGP braneworld cosmology and its family of mod- 
ifications of the Friedmann equation, also with freezing 
behavior, and the Albrecht-Skordis or exponential times 
polynomial potential, whose history cannot be classified 
as purely thawing or freezing. 

The dynamics is conveniently represented by the equa- 
tion of state, or effective pressure to density ratio, w, and 
its variation w' = dw/dlna = w/H where a is the expan- 
sion or scale factor. The Hubble parameter, or expansion 
rate, H = a/ a. We work in units where SttG = 1. 



A. PNGB Model 

Protected from radiative corrections by a shift sym- 
metry, this model possesses technical naturalness and is 
characterized by a symmetry energy scale / . The po- 
tential reads 



V{^) = K [l-f cos((/.//)] 



(1) 



with setting the overall magnitude, hence related to 
the present dark energy density. The equation of state, 
and the dynamics in general, is governed by / and the 
initial field position (f>i. (It is convenient, as seen from 
the form of Eq. ([T|), to use 4>i/ f instead of (pi.) 
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One can scan over the parameter space of these three 
variables and examine the evolutionary behavior and vi- 
ability as a dark energy model. Figure [1] shows a selec- 
tion of trajectories in the w-w' plane. The time coor- 
dinate runs along these tracks, and can be thought of 
as the scale factor a or the dark energy density fraction 
of the total energy density, 0^,(0). As we change or 
riu, = 1 — flm, where rim is the dimensionless present 
matter density, different points along a track for given / 
and 4>i/ f correspond to the present. In fact, for some pa- 
rameter values the dark energy never dominates and the 
density is restricted to Qw{o-) < ^w,ma.x < 1- One can 
show that for fixed then O^j^max cx; so models 

with symmetry energy scales much less than the Planck 
energy, f <^1, tend not to be viable. 

Figure [T] shows a wide selection of trajectories that 
reach — 0.72 at the present. They fan out across the 
phase space, including ones that lie outside the conven- 
tional thawing region 3(1 + w) > w' > 1 + w (although 
these start along w' — 3{1 + w) at early times) Q). The 
exceptions have / <C 1, and are not generic in that for 
/ ^ 1 we must fine tune ever more strictly the initial 
condition (t>i/f in order to achieve such a present den- 
sity. Figure [2] plots the allowed values of (j)i/f, which 
decrease rapidly, roughly as (0i//)max e^^/-^. For ex- 
ample, when / = 0.1, then the field must start exquisitely 
close to the top of the potential: rather than (f>i / f rang- 
ing freely over [0, tt], it is restricted to be less than 10^"^. 
For / = 0.05, this becomes (pi/ f < 10~^. Apart from un- 
naturalness, such values may run into physical problems 
such as a tachyonic instability [1, @, @ . 

In the future, the scalar field reaches the minimum of 
the potential and oscillates around it, giving an equation 
of state w — (matter-like) when averaged over many 
oscillations. We discuss this further in comparison with 
the next model. 



B. Linear Potential 

The linear potential tilts a flat potential, so the field 
rolls ~ although it is frozen by the large Hubble friction 
at early times. The potential is given by 

V{^)^V, + {^~<j>,)V' , (2) 

where V is the slope parameter, a constant 0, If 
the slope becomes too steep then the field never has time 
in its evolution to build up to appreciable energy den- 
sity before the kinetic energy becomes substantial and 
w > 0, causing the fractional energy density relative to 
the matter density to decrease with scale factor. The 
evolutionary tracks for this model fan out in the phase 
space within the thawing region (some examples for this 
and other models appear in i jllFl in Fig.[7|). 

Figure [3] shows the long time evolution of the PNGB 
vs. linear potential models, showing the similarity of the 
tracks at first, then the dramatic difference in the fate 
of the universe as the PNGB field oscillates, acting like 



f=0.2 



f = 0.5 
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FIG. 1: PNGB models fan out through phase space as 
their parameters vary (though still mostly within the thaw- 
ing region). At early times the models all start frozen at 
{w,w') = ( — 1,0) and thaw, with the scale factor increasing 
along each curve, although at different rates in each case. 
Here we end the tracks when = 0.72. 



matter in a time averaged sense, and the linear potential 
field shoots away, leading to deceleration and a cosmic 
doomsday collapse. 

We also consider the related family of polynomial po- 
tentials, V ~ with n — 2, 4:. These are also thawing 
models although their future behavior asymptotes to os- 
cillation about a zero potential minimum. Hence they do 
not runaway to negative potential and a rapid, dooms- 
day collapse. The equation of state during the oscillatory 
phase time averages to w — {n — 2)/{n + 2) 



C. SUGRA Model 

Tracking models have an early time attractor behavior 
that allows a large variety of initial conditions to give the 
same evolution in the matter dominated era, ameliorat- 
ing fine tuning of initial conditions fio'l . One example is 
the family of inverse power law potentials [Tl| . Includ- 
ing Planck scale corrections motivated by supergravity 
theory changes the potential to [l2| 

y(0) = Kr"e*'/2. (3) 

This has a local nonzero minimum, or cosmological con- 
stant. The equation of state behavior is governed by the 
power law index n. 
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FIG. 2: To achieve dark energy domination in the PNGB 
model before the field relaxes to its minimum, the initial field 
value 4>i must be small enough to give a long period of cos- 
mological constant-like, or frozen, behavior. For very steep 
potentials, i.e. low symmetry energy breaking scales /, the 
field must initially be extremely finely balanced near the top 
of the potential, with the curve showing the maximum 4>i/ f 
allowed to achieve fi„ > 0.72 at some point in the evolution. 



FIG. 3: Long term evolutions of PNGB and linear poten- 
tial models have distinct implications for the fate of the uni- 
verse, both different from the cosmological constant case. The 
PNGB field oscillates while the linear potential rolls to nega- 
tive infinity. Today the models shown have wo — —0.77, with 
w'o =0.47, 0.52 respectively. At a = 2, the PNGB model is 
still on the innermost track, with w — 0.11 (the curve end is 
at a = 6.6), while the linear potential is off the plot, on the 
way to collapse. 



Although the exponential factor has no effect on the 
attractor phase, it does permit a more rapid evolution 
after the field leaves that trajectory, moving the equation 
of state closer to w = —1. Since the dark energy has 
w = —2/(2 + n) while on the attractor trajectory, one 
requires n <C 1 in the inverse power law model to accord 
with observations; this is somewhat eased for the SUGRA 
model. Today the field properties can cover a wide swath 
within the freezing region 0.2w{l + 'w) < w' < 3w(l + iy). 



D. Braneworld Gravity Model 

Even dark energy theories that do not involve scalar 
fields can be viewed in terms of effective dynamics, where 
the equation of state is defined in terms of the Hubble 
parameter H(a) and its modified Friedmann equation: 



1 d\n6H^ 

Weft ^ ~l~ -— , 

6 a ma 



(4) 



where 6H^ — H^/Hq — ilma~^. One example involving 
very different physics from scalar fields is the extension of 
gravity theory through extra dimensions. This can lead 
to a modified Friedmann equation and effective equation 



of state (Ulil 



= p,n{a)/3 + {l-n„,)H^{H/HoY 

-1 



1 



(5) 
(6) 



where pm is the physical matter density and a is a param- 
eter depending on boundary conditions between our four 
dimensional universe and the higher dimensional bulk 
volume. The best motivated model in this family is DGP 
braneworld gravity, corresponding to a = 1 pLSl . Il6l | . At 
early times the effective potential looks like an inverse 
power law [l|, with index n — 2a/(2 — q), and so has 
tracking behavior. At late times the field rolls asymp- 
totically to a halt at a finite value of both the field and 
potential, acting as a cosmological constant. Indeed the 
trajectories lie within the freezing region. 



E. Albrecht-Skordis Model 

A scalar field potential with greater complexity is the 
Albrecht-Skordis p7| . or exponential with polynomial 
prefactor, potential, motivated by string theory. This 
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has the form 



(7) 



in the notation of [18], with a more compact but equiva- 
lent notation being 



V{^)^V,{l + Act>'')e-^'^, 



(8) 



where we shift to (j) — ip — (3, showing that only three 
parameters enter: 14, related to the dark energy density 
today, A, and A. Away from (j) ^ this behaves like 
an exponential potential, a classic tracker, so the initial 
conditions are not very important [l^. Near = 0, the 
potential has a false minimum, so a field rolling through 
this region can have complicated dynamics, and indeed 
be trapped and oscillate about the nonzero potential min- 
imum, eventually relaxing to a cosmological constant. 

Figure m show trajectories for different parameter val- 
ues, illustrating the wide variety of possible behaviors. In 
addition. Fig. O plots the equation of state w{a) so one 
has another view of the damped, oscillatory evolution. 
Note that while the field sees an exponential potential, 
away from the false minimum, it exhibits not only track- 
ing but tracing behavior - the dark energy equation of 
state is equal to the background, e.g. matter dominated, 
equation of state wi,. This means that the dark energy 
density is then a constant fraction of the background den- 
sity, given by flw,tiacc = 3(1-1- u;f,)/A^ So as not to 
violate primordial nucleosynthesis or cosmic microwave 
background constraints, this requires the contribution to 
be no more than a few percent. We show the dynam- 
ics for two cases, the first using the parameter values in 
p^ . corresponding to A = 3.4 and A — 106.7, which 
has an early dark energy fraction fie = 0.26 during mat- 
ter domination (fig = 0.35 during radiation domination), 
and the second using A = 10, keeping A the same, giving 
fie = 0.03 (0.04) during matter (radiation) domination, 
close to the upper limit allowed [H, HI] . 

For values of A allowed by nucleosynthesis and CMB 
limits, A > 10, the oscillations are absent or negligi- 
ble. One can show that the amplitude of the oscillations 
depends predominantly on the ratio X^/A (e.g. define 
(p = Xcf) and the potential only explicitly contains the 
parameter combination X^/A). If this combination ex- 
ceeds one, then there is no minimum but merely a slight 
local lessening of the exponential slope, and hence no os- 
cillations. The amplitude increases as /A approaches 
zero. However, since A > 10, small values of X^/A re- 
quire A > 1000, seemingly unnatural. Furthermore, the 
period of the oscillations is given by the effective mass 
and is inversely proportional to A for fixed X^ /A, and so 
for allowed A the oscillations will be negligible for z < 3. 
These behaviors are illustrated in Fig. [51 



F. Cross Comparison 

To compare the behaviors of different families, we plot 
selected representatives in Fig. [71 Varying the param- 
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FIG. 4: Albrecht-Skordis model acts like a tracer at early 
times, with a constant energy density fraction fig and w = 
in the matter dominated era, before oscillating around the 
nonzero minimum of the potential. We change the line thick- 
nesses (and colors) at z = 2, 1, and 0. The oscillations are 
invisible for more viable Qe ~ 0.03 case. 



eters within each model, as well as considering differ- 
ent models, spreads the evolution over regions of the w- 
w' phase space. Generally we see both similarities and 
distinctions between models and between families. One 
must also take into account the time coordinate along 
the curves, so that crossing of trajectories does not mean 
they have identical properties at any one moment. We 
plot the trajectories up to when the dark energy density 
is = 0.72. 

We could extend the curves into the future, as was 
done in Figs. [31 and [H The (/)^ potential will eventually 
settle dX w — 1/3, acting as radiation, as discussed in 
i jll Bl after oscillating around the minimum. Note that 
the SUGRA, DGP/iJ", and Albrecht-Skordis models all 
have nonzero minima, i.e. hidden cosmological constants, 
so they settle to w — —1. The SUGRA field does not os- 
cillate around the minimum because it approaches it with 
low kinetic energy, freezing to the cosmological constant 
state; the H" family has only an asymptotic minimum, 
also approached by freezing. 

As an alternative to showing the evolution of varied 
models at all times, we can take a slice at a particularly 
time, say when $7^, — 0.72, and construct phase space 
curves where the parameters of a potential vary along the 
curve. This can clear the illusion of overlap in behavior 
and provides an intermediate step toward the calibration 
in the next section. 




FIG. 5: Equation of state w{a) shows an alternate view of 
the evolution in Fig. [l] Note that for model parameters that 
do not violate early matter domination, the behavior relaxes 
swiftly to a cosmological constant, as shown by the solid, black 
curve. 



Figure [5] gives an example of this for the PNGB model, 
where the parameter running along the curve is the initial 
field position ipi/ f. That is, every point along any curve 
has riu, = 0.72 today, but corresponds to a different set of 
parameters for the potential and a different evolutionary 
behavior. This illustrates that different symmetry energy 
scales / define distinct paths to achieving a given dark 
energy density. (Of course not all of these are viable, with 
large values of along each curve corresponding to w 
far from —1, and small values of / suffering from the ex- 
treme fine tuning problem discussed in ijll AH However, 
by evaluating the equation of state and its time variation 
at a single time, we lose all dynamical information. In the 
next section we combine the advantages of the parameter 
scan with those of the evolutionary trajectories. 



III. STRETCHING DARK ENERGY 

To keep the dynamics central, we want to preserve in 
some way the temporal information, i.e. the field evolving 
from its high redshift state along a trajectory describing 
the equation of state and its time variation. However, we 
are free to rescale the time coordinate and define a time 
variation other than w' = dw/dlna. In particular, we 
can ask whether there is a global transformation that in 
some way calibrates the dark energy characteristics. We 
call this the evolutionary stretch factor. 



FIG. 6: The amplitude of oscillations in the equation of state 
is governed by X^/A and the period goes as (for fixed 
X^/A). The figure shows the equation of state history for 
A = 10 {Qe = 0.03 during matter domination) as solid curves, 
for /A = 0.1, 0.5, 0.9, from highest peak to lowest, and for 
A = 3.4 (fie ~ 0.26), with X^/A — 0.1, as a dashed curve. 
For appreciable oscillations X^ /A must approach zero, but for 
allowed (large) values of A any oscillations damp away for 
z<3. 



Stretching the time variation by different amounts at 
different times effectively introduces additional evolution 
beyond the scalar field behavior, so we consider a con- 
stant stretch factor, a simple renormalization. That is, 
we take w'{a) — > w'{a)/ai,. Now, since realistic obser- 
vations cannot map out the detail of the equation of 
state function, we seek to condense the information on 
the evolution to a set that is robustly constrained by 
data. Overcompression loses important physical prop- 
erties while undercompression leads to uninformatively 
large uncertainties. In the next section we will test the 
full stretch prescription to ensure that neither case oc- 
curs. To begin with, consider evaluating our new time 
variation quantity at a particular scale factor; further- 
more, to keep the number of parameters in the stretch 
prescription to a minimum, we choose this scale factor 
to be the same as the stretch factor a*. That is, the 
procedure can be viewed illustratively as 



For evaluating the value of the equation of state func- 
tion itself, w{a), we also avoid choosing an arbitrary scale 
factor. This leaves us with two choices: either a* or the 
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FIG. 7: Representative models considered in this section are 
plotted for various parameter values in the w-w' phase space. 
Solid, black curves are PNGB, short-dashed, red curves are 
for the linear potential, dotted, blue curves are for 4>^, long- 
dashed, black curves for the braneworld model (a = 1 DGP 
and a = 0.5), and dot-dashed, red curves for SUGRA. 



FIG. 8: These curves for the PNGB model correspond to a 
scan over the potential parameter space to find those values 
where fi^, — 0.72. Each curve is for a different energy scale /, 
with the parameter ^i// running along each curve, from zero 
at «; = — 1 to a maximum possible value shown in brackets. 



present epoch, a = 1. If we choose a*, then this proce- 
dure merely chooses a single point along the evolutionary 
trajectory, losing much of the global information. Thus 
we adopt wq = w{a = 1) and examine the dark energy 
characteristics in the plane of the two parameters, wq 
and w'{ai,)/ai, (= dw/da{a^))^ to see if there is indeed a 
normalizing relation for the time evolution. 

Figure [S] shows clearly that this prescription calibrates 
the evolution of the PNGB model. Instead of the fan of 
trajectories spreading through the w-w' phase space, as 
in Fig. [HI we now have a tightly calibrated, one parameter 
relation in the wq vs. w'{ai,)/ai, plane. Despite scanning 
over the model space of / and (pi/ f, this stripe is nar- 
row and well defined. Points within the stripe represent 
individual realizations of the PNGB model with choices 
of the symmetry energy scale ranging over the physically 
reasonable range / S [0.2, 5] and initial field position cov- 
ering from to the maximum value that allows Sl^, « 0.7. 

This tight calibration spreads little if we vary the 
present dark energy density as well as the potential pa- 
rameters themselves. Allowing VL^ today to range over 
0.69-0.75 gives the slightly wider, hghtly shaded region. 

Calibration succeeds for the other dark energy models 
considered as well, covering a wide range of physical ori- 
gins. Indeed, all the thawing models are closely related, 
nearly forming a single family under the calibration. The 
similarities extend to defining a single stretch parameter 



a* = 0.8 for the entire thawing class. Freezing fields 
also can be calibrated, with a uniform stretch parameter 

— 0.85, though the families stay more distinct within 
the freezing class. Figure fTOl shows the tight relations of 
the different dark energy models, in strong contrast with 
the "fan" nature of Fig. [T] 

From the form of the stretch calibrated time variation, 
w'{ai,)/ai,, we can recognize this as nearly identical to 
Wq, the dark energy variable in standard use, defined by 
[24| as Wa — — w'(a* — 0.5)/0.5 to fit the equation of 
state function by w{a) = wo + w'^\l — a). The super- 
script w indicates that the value of was chosen to fit 
w{a). Here, however, we defined the equivalent of Wa to 
calibrate dark energy families. This resulted in a* — 0.8 
for the thawing class and a* = 0.85 for the freezing class. 
An interesting further implication is that the "new" form 

wia)^wo-^^^^^{l-a)^wo+wlfHl^a) (10) 
a* 

has excellent accuracy when fitting the observables of 
distance and Hubble parameter, as we discuss next. 

IV. OBSERVING DARK ENERGY 

While the form (fTO|) was just shown useful in inter- 
pretation of dark energy theory, we should also investi- 
gate its utility for interpreting dark energy observations. 
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FIG. 9: Defining a new time variation variable Wa from 
w' calibrates the PNGB model into a tight locus; compare 
the spread in Fig. [S] Solid lines are for fixed / parameter, 
the shading shows the range of behaviors for / G [0.2,5]. 
The lighter shading shows the effect of also scanning over 
n,n = 0.25-0.31. 



Three related, but slightly different, physical bases exist 
for using the form w{a) = wq + Wa{l — a) to character- 
ize dark energy: this can be interpreted as 1) a fitting 
formula to the equation of state, 2) a calibration relation 
for families of dynamics in the w-w' plane, or 3) a fitting 
formula for observables such as distances and the Hub- 
ble expansion rate. The last two in particular are closely 
related and give similar results; indeed, when models do 
not deviate greatly from cosmological constant behavior 
the two approaches are almost identical. 

We now explore the accuracy of the form (fTU|) in fit- 
ting the exact distance-redshift and Hubble parameter- 
redshift relations for the diverse dark energy models dis- 
cussed in []IT1 

Table |T] summarizes the accuracies on d and H for a 
diverse range of models. These are generally good to the 
10"'^ level. Models closer to A would have better fits 
than shown here; models further from A are not favored 
by current data. For simplicity we henceforth denote the 
calibrated fit parameter simply as Wa- 

We could push the accuracy even further by mini- 
mizing the deviation not globally, over the entire range 
a G [0,1], but over a particular epoch, say a € [0.5,1]. 
However, we retain the global fit in general. Also, we 
have not taken advantage of the degree of freedom of 
wq, which could improve the fits. We emphasize that 
the stretch factor is a function of the dark energy physics 



FIG. 10: In terms of the calibrated dark energy parameters 
Wo and models and families lie in tightly homogeneous 

regions, in comparison to Fig. [T] showing the same models 
before calibration. We here vary over all parameters in the 
potentials. Shading shows the effect of scanning over ±0.03 in 
flm (we omit the shading for and linear potential models 
to minimize confusion; the width would be about half that 
shown for PNGB). Thawing models, despite their differences 
in w-w', are nearly identical once calibrated. Distinctions 
from freezing models, and between freezing models, become 
highlighted with calibration. 



and not dependent on the experiment, priors, etc. (in dis- 
tinction from a pivot redshift or pivot equation of state 
value). 

Note that the results from this prescription also an- 
swer the important question of whether the calibration 
procedure preserves the information faithfully to the pre- 
cision level of the data, or over- or under-compresses the 
model characteristics. A one parameter approach such 
as a constant value of w would have errors of order 1-2% 
in distance and up to 3% in Hubble parameter for the 
models we considered. This is insufficient for forthcom- 
ing observations. Conversely, since the two calibrated 
parameters of wq and Wa map the observables to better 
accuracy than expected from next generation data, these 
two parameters suffice and the data precision does not 
call for further equation of state parameters. 

This is not to say that some models could not exist 
where a third parameter carries information, but such 
models may not be generic or natural; the wide range of 
models considered here has no use for one. If we reach 
the stage of probing the cosmic expansion history below 
the 10~^ precision level, we should revisit the question of 
a further calibration parameter. 
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Model 


Sd/d 


5H/H 


PNGB {wo = -0.85) 


0.05% 


0.1% 


PNGB {wo = -0.75) 


0.1% 


0.2% 


Linear Pot. (too = —0.85) 


0.05% 


0.1% 


Linear Pot. {wq = —0.75) 


0.1% 


0.3% 


0* (wo = -0.85) 


0.01% 


0.04% 


(1)^ {wo = -0.75) 


0.02% 


0.06% 


Braneworld {wq = —0.78) 


0.03% 


0.07% 


SUGRA {n = 2) 


0.1% 


0.3% 


SUGRA (n = 11) 


0.1% 


0.3% 


Albrecht-Skordis {Q^ = 0.03) 


0.01% 


0.02% 


Albrecht-Skordis (fie = 0.26) 


0.1% 


0.4% 



TABLE I: Accuracy of Wo-Wa in fitting the exact distances 
and LIubble parameters for various dark energy models. These 
numbers represent global fits over all redshifts (except for the 
last three cases, where the fit covers z = 0-3, due to early dark 
energy: see 33). Better fits can be found over finite redshift 
ranges. 



Finally, this prescription is meant to help us find our 
way through the dark forest [2^ of models of cosmic ac- 
celeration, making accurate, more or less model indepen- 
dent assessments. Once precision data exist, they should 
be analyzed for every model of interest and within every 
applicable fit technique, parametric and nonparametric. 
We have seen that until we reach that point w^-Wa serves 
as a robust indicator and guide for predicting and com- 
paring cosmological probe information. 

V. FIGURES OF MERIT 

The accuracy of the wq and Wa form, defined in the 
manner discussed here, for characterization of observable 
properties of dark energy is at a level of order 10"'^, suf- 
ficient for next generation experiments. The calibration 
into tight families of equation of state properties, as seen 
in Fig. [3 suggests that not all combinations of wq and 
Wa are of equal insight. For example, one might distin- 
guish models in the thawing class from the cosmological 
constant and from each other by constraining the combi- 
nation varying exactly along the calibrated curve. 

Since this curve is nearly straight, we can characterize 
it by slope mj, and define a new variable 

w\ ='Wq + mtWa , (11) 

where the derivative with respect to this parameter runs 
parallel to the calibrated curve. Hence, determining 
localizes the behavior and distinguishes the specific dark 
energy characteristics. The narrowness of the calibrated 
region means that it is not so useful within thawing mod- 
els to constrain the direction perpendicular to the curve. 

One can define a similar variable for the freezing class, 
although here the families are more spread out, so the 



slope is more of an average than a well defined value, 

wjl = Wo + rufWa ■ (12) 

Values of mt — —1.75 and mj — 3.5 are reasonable 
choices. Note that the combinations wj| and w|| are not 
orthogonal, so the variable defined for each class does 
have utility in constraining the other class as well. For 
example, along the PNGB curve of Fig. [TO] the param- 
eter runs from —1 to 0, while goes from —1 to 
-2.2; along the SUGRA curve wf| = -1 to -1.7 while 
= — 1 to -1-0.9. This shows that each parameter, while 
optimized for a given physics question, does carry infor- 
mation on the other class. 

Thus, knowledge of either parameter w| or wj| an- 
swers the key questions of distinction from a cosmological 
constant, distinction between models, and to an extent 
distinction between classes. Constraining both parame- 
ters tightens the distinguishing ability, especially between 
classes, and provides a crucial crosscheck of the frame- 
work. 

It does not seem natural or effective to combine the 
uncertainties in estimating these variables from obser- 
vations into a single number, e.g. cr(i(;|) x a{w^^^), since 
they represent very different physics. Moreover, further 
investigation is needed into the optimum values for my, 
TOt and other issues before defining ultimate figures of 
merit, if this is even possible. However, the tightness 
of the calibration does imply that some combinations of 
Wq and Wa will provide insight into the nature of dark 
energy. Therefore, knowledge of the uncertainties cr(wo) 
and a{wa) and their covariance are the main ingredients 
for a variety of future figures of merit that might be de- 
veloped. 

Finally, we note that the accuracy of the w^-Wa form 
does start to degrade to the 10~^ level as dark energy 
becomes increasingly important in the early universe 
around z ^ 10'^, upsetting standard matter domination. 
See, for example, the last three models in Table U where 
the dark energy equations of state at recombination are 
w « —0.15, 0, 0, respectively. It could be useful to treat 
such early dark energy models as a separate class, and 
include constraint on the dark energy density lie at re- 
combination (which can best be done through growth 
probes) as another desideratum for a dark energy science 
program. 

VI. CONCLUSIONS 

Having investigated a diverse group of dark energy 
models to explain the acceleration of the cosmic expan- 
sion, we find a homogeneous "stretch" relation that cal- 
ibrates the time variation behavior into tight families. 
This stretch factor is closely related to the standard time 
variation measure Wa , and we verify that the equation of 
state form w{a) = wq + Wa{^ — a), with Wa now treated 
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as a fit parameter to observables, delivers fractional ac- 
curacy at the 10""^ level. 

Such accuracy is sufficient for next generation data and 
the WQ-Wa form can be viewed as an appropriate com- 
pression of the expansion history information that can 
be extracted from such observations. That is, this form 
neither overcompresses (loses important information) nor 
undercompresses (lacks additional leverage). This indi- 
cates there is no need nor generic benefit for going to a 
third parameter. Note that [1^ saw similar compression 
and tight relations within a principal component analysis 
relying on many modes. 

To gain insight into the nature of dark energy, par- 
ticular combinations of wo-Wa may have enhanced lever- 
age and hence merit, separating the cosmological con- 
stant from the thawing class, each from the freezing class, 
and possibly zeroing in on specific models within a class. 



The calibration, and its robustness and accuracy in ac- 
counting for the observable relations, offers a well-defined 
method for assessing the next generation dark energy 
science program. Interpretation of those observations 
should offer promising insights into the physics of the 
accelerating universe. 
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